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Abstract 

We use the chiral quark model to construct the complete 0{p^) weak A.S = 1 chiral 
lagrangian via the bosonization of the ten relevant operators of the effective quark 
lagrangian. The chiral coefficients are given in terms of /tt, the quark and gluon 
condensates and the scale-dependent NLO Wilson coefficients of the corresponding 
operators; in addition, they depend on the constituent quark mass M, a parameter 
characteristic of the model. All contributions of order N"^ as well as and agNc 
are included. The 75-scheme dependence of the chiral coefficients, computed via 
dimensional regularization, and the Fierz transformation properties of the operator 
basis are discussed in detail. We apply our results to the evaluation of the hadronic 
matrix elements for the decays — > 27r, consistently including the renormalization 
induced by the meson loops. The effect of this renormalization is sizable and intro- 
duces a long-distance scale dependence that matches in the physical amplitudes the 
short-distance scale dependence of the Wilson coefficients. 
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1 Introduction 



Chiral perturbation theory provides a faithful representation of the hadronic 
sector of the standard model at low energies. The form of this effective field theory 
is determined by SUl{3) x SUr{3) chiral invariance and its spontaneous breaking 
that, together with Lorentz invariance, dictate all possible terms. 

A particular example of chiral perturbation theory is the weak chiral lagrangian 
responsible for the AS = 1 transitions. This lagrangian controls most of the physics 
of kaon decays and in particular the AI = 1/2 selection rule and the CP- violating 
parameter e'/e, the determination of which is at the origin of this work. To the 
order 0(|?^), it is given in section 4. Such a lagrangian, by including also the effects 
of the isospin-symmetry violating electromagnetic interactions, is more general than 
that usually discussed in the literature. 

The coefficient in front of each term contains the short-distance information. 
Since we do not know how to connect directly the quark and gluon degrees of freedom 
of QCD to the hadronic states, these coefficients have to be determined either from 
a comparison to the experiments or by means of some phenomenological model. It 
is the purpose of this work to provide a theoretical determination of them by means 
of the chiral quark model (xQM) approach ^, ^ which is briefly described in 
section 3. 

In the xQM, and similar models, the interaction among mesons proceeds only by 
means of quark loops (that explains why these models are often called quark-loop 
models). Their theoretical justification can be found, in terms of a more fundamen- 
tal picture, as the mean field approximation of an extended Nambu-Jona-Lasinio 
(ENJL) model Q effectively mimicking QCD at intermediate energies. 

Starting with the short-distance effective AS = 1 lagrangian (see section 2) 
written in terms of four-quark operators, the xQM- allows us to compute the contri- 
bution of each of these operators to the corresponding coefficient of the low-energy 
chiral lagrangian. In the process, all the arbitrary coefficients are parameterized in 
terms of the input variables. These are the Wilson coefficients — that we take to 
the next-to-leading order (NLO) — of the corresponding operators, the pion decay 
constant /tt, the quark and gluon condensates and the constituent quark mass M, 
a parameter characteristic of the model. 

This determination requires some caution because of the presence of divergent 
contributions. Regularization and renormalization prescriptions must be given. In 
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order to be consistent with the regularization of the anomalous dimensions of the 
short-distance analysis, we use dimensional regularization in the modified minimal 
subtraction scheme. In this framework, we must deal with the problems of the 
75-scheme dependence and of whether Fierz transformations can be applied to the 
original basis of chiral operators. These problems are discussed in section 5, where 
we show that the t'Hooft-Veltman (HV) scheme generates some fake chiral anomalies 
that must be subtracted. 

We show that operators related by a Fierz transformation lead to equivalent 
contributions to the chiral lagrangian coefficients in the HV scheme but not in the 
naive dimensional regularization (NDR) one. For the latter, Fierz transformations 
with respect to the basis used for the calculation of the Wilson coefficients must be 
avoided. As an outcome, a new 75-scheme dependence arises in the coefficients of 
the chiral lagrangian in addition to that, already present from the beginning, of the 
NLO Wilson coefficients. The matching between the two, which depends of the input 
parameters, and the resulting 75-scheme independence, can be used to restrict the 
values of M, the only free parameter in the model. While such a procedure is only 
mentioned in passing here, it has been discussed in |5| and will be used extensively 
in future work. 

In section 6, we write explicitly all the chiral lagrangian coefficients determined 
in the HV and NDR schemes. Both schemes are consistent for our lagrangian that 
contains no anomaly. We have included contributions of order 0{N^), 0{Nc) as 
well as 0{asNc)- The latter represents non-perturbative corrections induced by the 
gluon condensate; their relevance was first advocated in [^] where the effect on the 
operators Qi and Q2 was computed. 

In section 7, we discuss the numerical values of the input parameters in terms of 
which these matrix elements are given by the xQM^- The determination of these input 
values is delicate and it requires some discussion because of the various estimates 
presented in the literature. 

For future applications, we specialize our results to the decays 27r and 

give, in section 8, all the relevant hadronic matrix elements to order 0{asNc). Such 
an order includes the leading (that is 0{N^)) as well as the next-to- leading order in 
the 1/Nc expansion and the leading corrections due to the gluon condensate. The 
reader who wishes to skip the technical details of how the result has been obtained 
is advised to go directly to this section. 
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The renormalization of the amphtudes induced by the meson loops (a correc- 
tion of order 0{Nc)) is sizable and it is given in Section 8.1. Here our approach 
differs with respect to conventional chiral perturbation theory in the way the scale 
independence of the physical amplitudes is obtained. Whereas in the conventional 
approach the scale dependence of the meson-loop renormalization is compensated 
(by construction) by that of the counterterms in the O(p^) weak chiral lagrangian, in 
our approach these counterterms are taken to be scale independent at the tree-level. 
The chiral loops provide the scale dependence of the hadronic matrix elements that 
is absent in the xQM itself. It is through the matching between this long-distance 
scale dependence and the short-distance one contained in the Wilson coefficients 
that the scale independence of the physical amplitudes should finally be achieved. 

A similar point of view was first advocated and explored in ref. |^] (and sub- 
sequently in [^) in the framework of a cut-off regularization for the meson- loop 
renormalization (while using dimensional regularization for the Wilson coefficients). 
Our approach, based on dimensional regularization of the meson loops, leads to 
different results, as discussed in section 8. 

In the appendix, the interested reader can find the relevant Feynman rules and 
other formula useful in the computation. We also briefly summarize the HV and 
NDR 75-schemes and give a table of the numerical value of all the input parameters. 

2 The AS* = 1 Effective Quark Lagrangian 

In our approach, the weak chiral lagrangian is determined by the ten operators 
defining the AS = 1 effective lagrangian at scales /x < nic |l^, |ll|: 

^A5=i = -^V^dV:,J2[z,ifi) + Tyi{f,)]Q,{f,) . (2.1) 
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The Qi are four-quark operators obtained by integrating out in the standard model 
the vector bosons and the heavy quarks t, b and c. A convenient basis is the following: 

Ql = (San^)v__A (%^a)v-A ' 
Q2 = (su)v-A M)v-A ' 
Q3,5 = (sc^)v-A Eg (^g)vTA ' 

Q7,9 = I (sd)y„A Eg eg (gg)v±A ' 

Qs.lO = i (Sad/3)v_A 5I]g eg(g/3ga)v±A , 

where a, j3 denote color indices = l,...,A'c) and eg are the quark charges 

(erf = Cs = —1/3 and e^ = 2/3). Color indices for the color singlet operators are 
omitted. The subscripts iV ±A) refer to the chiral projections 7^(lib75). We recall 
that Qi^2 stand for the Ty-induced current-current operators, Qs-e for the QCD 
penguin operators and Qt-iq for the electroweak penguins (and boxes). 

The functions Zi{^) and yj(/i) are the Wilson coefficients evaluated at the scale 
/i and Vij the Kobayashi-Maskawa (KM) matrix elements: r = —VtdV^g/VudVus- 
Following the standard parameterization of the KM matrix, we have that the Zi{fj,) 
control the CP-conserving part of the amplitudes while the yj(ju) the CP-violating 
one. 

The numerical value of these Wilson coefficients depends on a^, the recent de- 
termination of which from NLO calculations on data taken at LEP and SLC [^] 
gives 

a,(mz) = 0.119 ±0.006, (2.3) 
which corresponds to 

AgJ^P = 350 ± 100 MeV , (2.4) 

the range of values we will use for numerical estimates. 

Even though not all of these operators are independent, the basis in eq. (|2.2| ) 
has become the standard one. As a matter of fact, it is the only one in which the 
Wilson coefficients have been estimated to the NLO order |jll|, p^ ]. 

Two other possible operators, Qn and Q12, representing the dipole contribution 
of the penguin operators, only contribute to the next order in chiral perturbation 
and turn out to be negligible On the other hand, some of the other ten operators 



in eq. ( |2.2[ ) may receive sizable NLO corrections [14, 15, 



4 



3 The Model 



In order to evaluate the bosonization of the quark operators (|2.2| ) we exploit the 
xQM approach which provides an effective link between QCD and chiral perturba- 
tion theory. 

The xQM can be thought of as the mean field approximation to the ENJL model 
of low-energy QCD. A detailed discussion of the ENJL model and its relationship 
with QCD — as well as with the xQM^ — can be found, for instance, in ref. Q. 

In the xQM, the light (constituent) quarks are coupled to the Goldstone mesons 
by the term 

^xQM = -M (qji T^qi + u ^Ur) , (3-1) 
where = (n, d, s) is the quark flavor triplet, and the 3x3 matrix 

E = exp (^y n(x)) (3.2) 

contains the pseudoscalar octet n(x) = J2a A"7r"(x)/2, (a = 1, ...,8). The scale / is 
identified with the pion decay constant f^^ (and equal to fx before chiral loops are 
introduced). 

Under the action of the generators Vr and Vl of the chiral group SUr{3) x 
^[/^(S), S transforms linearly 

S' = Vr^vI ■ (3.3) 

with the quark fields transforming as 

Ql = QL and Qr = VRqR, (3.4) 

and accordingly for the conjugated fields. 

The lagrangian (|3.1|) is invariant under the chiral group. If we consider the field 



[/ = MS (3.5) 

as the exponential representation of the nonlinear sigma model, where M is invariant 
under SU{3)v and S is given by eq. (|3.2|) , the xQM can be understood as the mean 
field approximation in which M = MI. The quantum excitations associated to the 
full M are thus decoupled in a manner that preserves chiral symmetry. The quantity 
M, which is characteristic of the model, can be interpreted as the constituent quark 
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mass that is generated, together with the quark condensate, in the spontaneous 
breaking of chiral symmetry. It is given by 

M = mg- g—j- , (3.6) 

where 5 is a compHcated function of the modehng of the solution of the spontaneous 
breaking of chiral symmetry in QCD and 




A^ = 2^W— «0.82 GeV. (3.7) 



In our approach we have taken M as a free parameter to be eventually constrained 
by the 75-scheme independence of physical amplitudes (see the final comment of 
section 8). 

Fields other than the Goldstone bosons (like the p, for example) can in principle 
be included by going beyond the mean field approximation. They give well-defined 
and computable corrections [^. 

The gluon degrees of freedom of QCD are considered as integrated out down to 
the chiral breaking scale A-^, here acting as an infrared cut-off. The effect of the 
remaining soft gluons are assumed to be well-represented by the various condensates, 
the leading contribution coming from 

i—GG) . (3.8) 
vr 

The constituent quarks are taken to be propagating in the fixed background of such 
soft gluons. 

The xQM was first discussed in 0, |2|. New life was breathed into it in a series 
of more recent works [^]. 

We opted for the somewhat more restrictive definition suggested in |jl6| (and 
there referred to as the QCD effective action model) in which the meson degrees of 
freedom do not propagate in the original lagrangian. 



We have used the xQM of I'ef. [16| because: 



It is sufficiently simple without being trivial. It reproduces important features 
we know from the experiments like, for instance, the values of the O(p^) strong 
chiral lagrangian coefficients -Li's; 

It is sensitive to the 75-scheme dependence of dimensional regularization; 
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• By freezing the meson degrees of freedom, it allows us to deal with them in a 
separate step of the computation (see section 8). 

The model interpolates between the chiral breaking scale and M. The con- 
stituent quarks are the only dynamical degrees of freedom present within this range. 
The Goldstone bosons and the soft QCD gluons are taken in our approach as external 
fields. We neglect heavier scalar, vector and axial meson multiplets. 

A kinetic term for the mesons, as well as the entire chiral lagrangian is generated 
and determined by integrating out the constituent quark degrees of freedom of the 
model. The AS" = 1 weak chiral lagrangian thus becomes the effective theory of 
the xQM below the scale M. In the process, the many coefficients of the chiral 
lagrangian are parameterized — to the order 0{asNc) in our computation — in terms 
of the input variables f-,^, the quark and gluon condensates and M, the only free 
parameter of the model. 

In the strong sector, where the O(p^) coefhcients Li are experimentally known, 
these can be used to compare the prediction of the model. These coefficient have 
been computed in [|^. To the leading order in Nc, the Li (except for L5 and Lg) are 
purely geometrical factors that cannot be compared directly with the experimental 
values which have an explicit scale dependence. Only combinations of the same that 
have vanishing anomalous dimension can be compared. In this case the result is 
quite encouraging. See, for instance, L3, Li — L2/2 and Lg + Liq. A more stringent 
comparison can only be made by including the heavier multiplets that are known to 
give rather large contributions to these coefficients. A recent computation shows an 
impressive agreement Q]. 

In the ENJL model also the other parameters of the strong sector — like, for in- 
stance, M and /jr — are in principle computable quantities themselves. As attractive 
as this line of investigation might be, we do not pursue it here. In our approach, 
the strong sector parameters are only taken as input values, determined from the 
experiments, in order to predict the physics of the weak sector. 

We would like to emphasize here the novelty of our approach with respect to 
the usual treatment of the weak chiral lagrangian. This is best understood by 
comparison with the strong sector where the xQM uses the value of as an input. 
In this case the coefficient of the chiral lagrangian is scale independent and there 
is no matching to any independently determined Wilson coefficient. Because of 
this, the strong sector follows the usual treatment and in particular the higher-order 
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counterterms must have a scale dependence in order to cancel that induced by the 
chiral loops. This is done by construction. The weak sector behaves differently. 
The xQM determines the coefficients as a function of the Wilson coefficients so 
that they have an explicit scale dependence that should then be matched to that 
independently induced by the chiral loops. Whether this happens represents a crucial 
test of our approach. In this case the O(p^) chiral counterterms need not have any 
scale dependence at the tree level. 

As it is the case for the counterterms Li in the strong sector, the effect of the 
heavier multiplets on the xQM approximation may be important. In particular, it is 
known Q that the vector multiplets mix with the Goldstone bosons thus modifying 
the axial coupling qa between constituent quarks and mesons that in the xQM is 
taken to be 1. Their inclusion is postponed to future work. 

3.1 BOSONIZATION 

A simple procedure allows us to find the bosonization of the four-quark oper- 
ators in eq. ( p.2[ ). First, we reproduce their SU{?>) fiavor structure by using the 
appropriate combinations of Gell-Mann matrices acting on the quark fiavor triplets 
q. For instance, left-handed current operators (like Qi and Q2) are written as 

qiKl^qi qL>?gl^.qL ■ (3.9) 

where AJ^ (m, n = 1,2,3) are the appropriate flavor projectors. Then we rotate the 
quark multiplets as follows: 

qL = i'^Qh , qi = Qhi 

qR = CQr, qR = QRi\ (3.10) 

where = S and (C^)^ = S^. The Ql,r are the constituent quark flelds. 

The rotation ( |3.1[1| ) is such as to transform the lagrangian of eq. ( |3.1D into a pure 
mass term for the constituent quarks, 

-M^Qj^Ql + QlQr) , (3.11) 

while the quark coupling to the Goldstone bosons is transferred to the kinetic term 
in the QCD lagrangian. The same rotation, applied to the operators in eq. (|3.9| ), 
together with their Fierzed transformed expressions, yields 

QLi\nl^£}QL QLiKl^.£}QL and (n ^ (?) . (3.12) 
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Finally, the rotated quark fields are integrated out after having attached the 
axial fields 

A, = -y{D,J:^y = y^{D,Y:y (3.13) 

to the quark loops in all possible manners that are consistent with Lorentz invariance. 
In the case of eq. ( |3.12 ), the insertion of no axial fields (constant term) gives no 
contribution since 

Tr (CA-e^ eA^C^) = Tr (A^A^) = , (3.14) 
unless {p,q) = {n,m). The first non-vanishing contributions are proportional to 

Tr(^A™e^^^CA^e^^M) ^ Tr(A™StZ),,SA^StZ)^s) 
Tr {^X::^UX^,^^ A^An oc Tr (a™A^ Sti?^S StZ)^s) (3.15) 
Tr(CA-Ct^/^^^^A^^t) ^ Tr(A^A-StZ)^SStZ),,s) 

which provide, together with the (n q) expressions, the O(p^) chiral representation 
of the quark operator. The single trace of eq. ( 0.15 ) can be written as the product 
of two traces (see appendix A. 4); in this form, the independent terms are more 
easily recognized and it is easier to verify that the bosonized forms satisfy the CPS 
invariance |17] of the original quark operator. 

The bosonization of operators involving the right-handed currents proceeds, 
along similar lines. Starting from the rotated operator 

Ql^XuI^C^Ql QrC^X^^^,,^Qr, (3.16) 
or its Fierzed expression 

QlCX^^Qb. Qr^^KC^Ql , (3.17) 
we obtain after the identification of equivalent terms 

Tr (a;^ St XP S 

Tr (A™ Z)aSt XP D^'S , , 

> / / \ (3.18) 

Tr (a™ SD^St XP StD^^Sj , 

Tr f A^f St XP D^T.D'^T.^ S 

Not all of these bosonizations are actually present for each operator. For instance, 
in the case of the gluonic penguins, the sum over the quark flavors together with 
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unitarity make all but one of these terms vanish. It is only for some of the electroweak 
penguins that all contributions are actually there. 

The contributions arising from the last term in eq. ( p.l8| ) are usually not included 
in the literature (more on that in section 8.2). As a matter of fact, by considering 
the 0{p^) bosonization of the quark density 



qLQR 



X X 



(3.19) 



the neglect of the last term in eq. ( |3.18| ) corresponds to discarding the second one 
of the two quadratic terms ^ . In the first paper of ref . the authors claim that the 
constant C2 (ci) can be put to zero by a nonlinear transformation that preserves the 
unitarity of S. We have verified that such a transformation has no consequences in 
the case of the gluon Penguins, where all 0{p'^) terms lead to the same bosoniza- 
tion (proportional to the combination ci + C2), while it cannot be applied in the 
electroweak sector, where non-equivalent bosonizations are generated by the two 
independent terms (see section 8.2). 

4 The Weak Chiral Lagrangian 

The strong chiral lagrangian is completely fixed to the leading order in momenta 
by symmetry requirements and the Goldstone boson's decay amplitudes: 

e2 , . f2 



St'rU = ^ Tr (D^SD^St) + Lbo Tr (mt) + SA^t) , (4.I 



r , , , 

'-'strong ^ y fJ- — y ■ 2 

where M = diag[m„, m^, m^] and Bq is defined by {qiQj) = —f^Bo6ij. To the 
next-to-leading order there are ten coefficients Li to be determined |l^ . 

For the weak AS = 1 flavor changing interactions, the systematic application of 
the procedure described in the previous section to the operators in eq. ( ^ ) leads, 
after some algebraic manipulations, to the following O(p^) chiral lagrangian: 



4'Li= GW(Q,,8)Tv(AistAls) 

+ GsiQs-w) TV (A^Z)^StZ)^'S 



third term C3 Df_iTiD^Ti^ E can also be included. The three terms are however linearly 



dependent and it is sufficient to keep only two of them as in eq. (3.1!:). The equivalence of any 
choice of two terms can be explicitly verified. To this end it is important to remember that the 
term proportional to C3 induces a wave function renormalization that must be included in the 
determination of the coefficients, see section 8.2. 
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+ Gii(Qi,2,9,io) Tr (AfstD^s) Tr (A^St^^S 
+ Gli(Qi,2,9,io) Tr (A^StD^s) Tr (AjstD^S 
+ GiH(Q7,8) Tr (a?Z)^s) Tr (A^D'^St) 
+ ^1^(^7,8) Tr (Aist^^s) Tr (AjSZ^^St) 
+ ^2^(^7,8) [Tr (a?s) Tr (A^Z^^StD^S St) 

+ Tr (AfD^^SD^^St S) Tr {^l^^)] , (4.2) 

where A*- are combinations of Gell-Mann SU{3) matrices defined by (A* );^ = 5u6jk 
and S is defined in eq. ( |3.2| ). The covariant derivatives in eq. ( |4.2| ) are taken with 
respect to the external gauge fields whenever they are present. 

The convenience of the non-standard form of eq. (|4.2|) will become clear in the 



following. In our approach the lagrangian (4.2) is the effective theory generated by 
integration of the three light quarks of the AS" = 1 quark lagrangian. Therefore 
the notation is such as to remind us the flavor and the chiral structure of the quark 
operators. In particular, Gs represents the (8^ x 1^) part of the interaction, as it is 
induced in QCD by the gluonic penguins, while the two terms proportional to G^^, 
and G^j^ are admixtures of the {27 l ^ ii?) ^^'^ the (8^, x 1^) part of the interaction, 
such as it is induced by left-handed current-current operators; the term proportional 
to G^'^^ is the constant (non-derivative) part arising in the isospin violating and 
{^j^^^^r) electroweak components, of which the terms proportional to G'^jf'^ represent 
the 0{p'^) momentum corrections. 

The terms proportional to Gg, G\j^ and G\j^ have been already studied in the 
literature |2|, ^, |l^ in the framework of chiral perturbation theory. Those propor- 
tional to and are usually separated further into their isospin components 
(see the appendix) as 

"2 
3 



C 



27 



927 



- Tr (a^S+D^s) Tr (AfS+D'^s) 



+ Tr (\It)D^t) Tr (A^StD^S 



(4.3) 



which is a (27^ x 1^), and 
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Tr 



Tr (a^S+D^s) Tr (AfS+D^S 
(AiS+D^s) Tr (A}stL>^s) 



(4.4) 
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which is a pure (8^ x ii?)- We prefer to keep the AS = 1 chiral Lagrangian in the form 
given in eq. (|4.2| ), which makes the bosonization of the various quark operators more 
transparent, and perform the isospin projections at the level of the matrix elements. 
Equations ( |4.3D ~(4.4) provide anyhow the comparison with other references. The 
chiral coefficients in the two approaches are related by 

1 



58(Ql,2) 
g27iQl,2) 



3G2i(Qi,2)-2Gl^(gi,2) 



GlLiQl,2)+GidQl,2) 



(4.5) 



Concerning the (8^, x 8^) part, the constant term was first considered in [^]; its 
momentum corrections are discussed here for the first time, as far as we know. 
In this paper we only need the weak chiral lagrangian to 0{p'^). The O(p^) weak 



lagrangian is very complicated with thirty-seven [14| coefficients to be determined 
only in the (8^ x 1^) sector, and many more in the others. Terms proportional to 
the current quark masses belong to such a higher-order lagrangian — except for a 
term 0{mq) for the operators Qj^s'- 

G(™)(Q7,8) [Tr (A^StA^STWts) + Tr (AjSA^St^WSt)] , (4.6) 

whose contribution is small compared to the leading constant term. 

The weight of some of the next-to- leading order corrections on the physical ampli- 
tudes has been estimated in various models |l^ as well as in the xQM They 
range from 10 to 30% of the leading contributions, a potentially large effect that we 
should bear in mind when estimating the inherent uncertainty of our computation. 



5 The Model at Work 

The seven coefficients Gi in ( |4.2D would have in general to be determined phe- 
nomenologically. In the xQM, they are calculable by integration of the constituent 
quarks. This section contains those details of the computation that require some 
explicit discussion, since a few subtle points arise in the regularization procedure. 



5.1 Regularization 

Some of the constituent-quark loops are divergent. We use dimensional regular- 
ization (with the definition d = 4 — 2e). The logarithmic and quadratic divergencies 
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of the loop integration are identified with, respectively, the pion decay constant /^r 
and the quark condensate (qq). We define the two quantities 

and then identify f^^^ = f-^ and (qq)^^^ = (qq) in the end of the computation. 
The latter are understood as the physical quantities, inclusive of gluon and mass 
corrections. In particular, a direct computation within the model shows that 

(gg)*-"^-' = (qq)^^'^ — \ (— GG) + (higher gluon condensates) , (5.3) 
12m vr 

where we neglect current quark mass corrections. Similarly for f^^\ We use 
eqs. ( |5.1| )-( |0| ) as a convenient bookkeeping device for these input parameters. 

There is a certain degree of arbitrariness in these definitions inasmuch as the two 
gamma functions go into each other up to a finite term by expanding around the 
pole. Yet, for all practical purposes, keeping the two types of divergencies separated 
is a simple and effective way of singling out the different contributions. 



Since the fl dependence in eqs. (5T) and (5^) is absorbed in the physical quan- 
tities fjr and (qq) |^, there is no scale dependence induced by the xQM in the chiral 
coefficients. 

The xQM can be thought as an effective QCD model interpolating between 
and M. In this sense, the most natural regularization scheme is a cut-off theory 
where no divergencies arise. In fact, one finds 

^ 47r2/ M2 ' ^""^^ 



Eqs. ( |5.4| ) and ( |5.5| ) are finite (albeit cut-off dependent) quantities to be identified 

with and (qq), respectively. Eq. ( |5.5| ) can be understood in the ENJL model as a 

solution of the Schwinger-Dyson equation with gap M. Notice that in this case M 

exhibits an intrinsic scale dependence and it vanishes at energies higher than A^^,. 

^The quark condensate has a perturbative scale dependence which originates in the short-distance 
computation; see the discussion of section 7 



13 




+ 







A 



B 



Figure 1: The quark constituent loops for an arbitrary operator. (A) is the unfactor- 
ized pattern, (B) the factorized one. The crossed circles represent operator and/or 
current insertions. 

Unfortunately, such a cut-off regularization of the hadronic matrix elements is 
not consistent with that of the Wilson coefficients, that are obtained by means 
of dimensional regularization. For this reason, we do not pursue this possibility 
further (see, however, refs. @, where a matching of cut-off regularization and NDR 
was attempted). 

The scale dependence — that must appear in the hadronic matrix elements in 
order to make physical quantities scale independent — is introduced in our approach 
by the meson loop corrections to the hadronic matrix elements, that we calculate in 
section 8 for the K tttt amplitudes. 

5.2 The Computation 

Within the model, any computation is easily performed by applying the Feynman 
rules given in the appendix. 

As shown in fig. 1, for each operator there are two possible configurations that 
must be estimated. Whether the unfactorized form of fig. 1(A) or the factorized 
form of fig. 1(B) is the leading one in the 1/Nc expansion depends on the color 
structure of the operator. For example, the unfactorized configuration is the leading 
one for the color "unsaturated" Qi, but the opposite is true for the color "saturated" 
operator Q2- 

The coefficients G^'^, for example, are most easily computed by analyzing the 
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A B 

Figure 2: The constituent quark loops coupled to the mesons for an arbitrary op- 
erator. (A) is the unfactorized pattern, (B) the factorized one with the insertion of 
the gluon condensate (in this case only the color octet components, proportional to 
T"", of the currents contribute). Meson and gluon lines are attached in all possible 
ways. 



two-meson diagrams of fig. 2. In the case of Qi, we may use the transition ir^ 
in order to fix uniquely G\i^{Qi) and —>■ 7r+ to fix G\]^{Qi) (see the Feynman 
rules of table 1 in appendix B). In the HV scheme and before gluon corrections they 
turn out to be 

GldQi) = -f\f^'^f/Nc, (5.6) 

which is the expected pattern in the expansion. 

The diagrams with three meson external lines, which are the relevant ones in 
the computation of the matrix elements of section 8, are then generated by means 
of the chiral lagrangian in eq. (|4.2|). 



5.3 Non-Perturbative Gluonic Corrections 

Non-perturbative gluonic corrections are introduced by propagating the quarks 
in the background of external soft gluons, as suggested in ref. Q. They provide a 
sizable correction of order 0{asNc) that we parameterize in terms of the leading 
gluon condensate {asGG/n) by introducing the quantity 

^(GG) = ^^^^. (5.7) 
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Since this correction is finite we can always compute it in the factorized form 
by means of a Fierz transformation (see eqs. ( A.ll| )-( |A.12| ) in the appendix). In 



this way it is easier to single out those contributions that are non-vanishing in the 
presence of the external gluons. The leading gluon condensate contribution is of 
0{1/Nc) and it is generated by diagrams of the type shown in fig. 2(B). 

Only those configurations in which one external gluon field is attached to each 
quark loop are genuine corrections to the matrix element, since those in which both 
gluons are attached to the same quark loop are included in the renormalization of 
either the quark condensate or /^r, for which we take the physical values. 

5.4 75-ScHEME Dependence 

In dealing with the problem of 75 in dimensional regularization, we have con- 
sidered both, the NDR scheme — anti-commuting 75 in d dimensions — and the HV 
scheme — commuting 75 in d 7^ 4 dimensions. Both procedures have been used to 
obtain a consistent set of NLO anomalous dimension matrices for the perturbative 



evolution of the Wilson coefficients 1 1 1 , ^ . 

Two questions related to the 75-scheme dependence need addressing before pro- 
ceeding with the actual computation. To these we now turn. 

5.4.1 Fake Anomalies in the HV scheme 

A first problem arises in considering the building blocks out of which the matrix 
elements in the factorized form of fig. 1(B) are made. These have an independent 
physical interpretation, as we shall see, and come in four kinds. Let us consider first 
the result of the NDR scheme. There are two densities, that we need to O(p^), for 
which we obtain 



(0|s757x|K+(fc))j,oR = iV2 



f 2M 



(5.8) 



and two currents, that we only need to 0{p): 

(0|s7^75«|i^+(fc))NDR = iv^A;^/^°^ (5.10) 
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(^+(p+)|^7^d|i^+(fc))NDR = + (5.11) 

where q = k — pj^ and P = k + while /^^^(O) = f^^^/f is identified with the 
vector form factor at zero momentum transfer q. 
In the NDR scheme we correctly find 

/^i^^O) = l and /-(0) = 0. (5.12) 

This result holds in the limit of unbroken SU{3) symmetry, where 

U = fK and /f"-"(0) = /r-°(0)=/±(0). (5.13) 

In this limit, eqs. (5.12) are in agreement with the experiments (deviations of /+(0) 
from unity are of order ||2^) and we do not find any anomalous result in the 
NDR scheme. 

As we consider next the HV scheme, the two densities are now given by 



{0\s-f5u\K+{k))^^ = (0|s75n|i^+(A:))NDR 

+ iV2 f 



12^ I 1 



A2 I 6M2 



(5.14) 



(^+(p+)|^d|K+(^))HV = (vr+(p+)|sd|/^+(A;))^j3j, + 24 — . (5.15) 



A2 



while for the current matrix elements we find 



(0|s7^75M|E:+(/fc))jjv = iV2kf'f'-^^ (5.16) 
{7r+{p+)\srd\K+{k))ny = - {q^)P^ + Uq^)q^ , (5.17) 



/HV(o) = i + 4_, (5.18) 



where 

X 

Eq. ( |5.16| ) is the same as eq. ( |5.1C1| ), so that f^^ is defined identically in the two 
renormalization schemes. On the other hand, the vector form factor /+(0) turns out 
to be different, as it happens for the density matrix elements. 

Even though {&) ^ 1 the vector current Ward identity is preserved also in 
the HV scheme. In fact, in the model the mesons propagate via quark loops, so that 
a simple calculation in the NDR case leads to 

G^\k)=k''ff'\^)- (5.19) 
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while in the HV we find 



24 

X 



(5.20) 



By looking at the term proportional to k'^ we see by inspection that the same shift 
in /+(0) is present in the propagator as well as in the vertex. Therefore the vector 
Ward identity 



1 



d^^d\7r+ip,))=G-'ip2)-G-'ipi) 



(5.21) 



holds in both schemes, and one might think that a redefinition of /^^(O) could solve 
all problems. Unfortunately, because /^^^ is the same in the two schemes, it is not 
possible to simply redefine f¥^{0) to be equal to 1, reabsorbing the HV shift in f^^. 



Moreover, eq. (5.20) shows another problem: a mass term is generated in the HV 
case, thus leading to explicit breaking of chiral symmetry f\. 

We therefore find that, in order to maintain the symmetries of the theory, we 
must subtract appropriate terms in the HV case. In particular, the mass term in 
eq. (5.20) must be taken away. This subtraction leads to an analogous subtraction in 
the building blocks of eqs. (|1D-(|1|), which become identical to the NDR results. 
Having fully subtracted the propagators, the Ward identity in ( p. 21 ) implies that 
also /^^(O) in eq. ( 5.17| ) becomes identical to the NDR result. 

The overall effect of the subtractions is to make the building blocks of the HV 
scheme identical to those of NDR, as given in eqs. ( ^ 



dSTlD. 

These subtractions can be implemented in the strong sector by the addition of 
appropriate terms in the xQM one-loop chiral lagrangian that becomes, in the HV 
case: 



C 



where 



(2) 

strong 



ai Tr 



ai 



a2 



a2 Tr 



X 



(5.22) 

(5.23) 
(5.24) 



^ There has been some discussion in the hterature about similar problems in the standard model. 
For instance, in ref. chiral anomalies in addition to the ABJ's were noticed. In ref. the 
need of subtractions was pointed out and some claims of ref. [pi| adjusted. 
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The lagrangian (|53^ ) suffices in giving the correct propagator to one loop. Notice 



that the term proportional to 02 breaks explicitly chiral invariance. 

Other terms, containing the appropriate flavor projectors, must be added in 
order to eliminate fakes anomalies arising in the weak sector. 

5.4.2 FiERz Transformation in the two schemes 

Another related problem is the following. By a Fierz transformation it is possible 
to bring all matrix elements (leading as well as next-to-leading order in 1/Nc) to 
a factorized form. As shown in the previous subsection, the factorized building 
blocks (after subtraction) are scheme independent. Therefore it would seem that 
there is no 75-scheme dependence at all in the hadronic matrix elements. Yet, 
this conclusion is not correct because, while in the HV scheme it is possible to 
apply a Fierz transformation on the quark operators without changing the matrix 
element, this is not possible in the NDR scheme where operators related by a Fierz 
transformation lead to different matrix elements. This fact is at the origin of the 
75-scheme dependence of the matrix elements. 

As a first example, let us consider the contribution of Qi to G^^. A direct 
evaluation of the unfactorized diagram in fig. 3(A) gives, in the NDR scheme, 

GUQi) = -fif^'^f{l-Q^^ ■ (5.25) 
In the HV case we obtain instead 

G\l{Qi) = -fif^'^f (5.26) 

We could have performed the same computation by first Fierz transforming the 
operator Qi that becomes 

Qi = {sd)v-A{uu)v-A- (5.27) 



In this case, by means of the building blocks in eqs. (|5T0|)-(|5Tl|) (or the subtracted 



eqs. ( |5.16| )-( |5T7 )), we would have obtained the result in eq. ( 5.26 ) in both the NDR 



and HV schemes. We explicitly see that the Fierz related operators Qi and Qi might 
lead in the NDR scheme to different contributions to the chiral coefficients. 

As a second, more complicated example, let us take the contribution of the 
operator Qq to Gg and evaluate its bosonization directly in the form given in eq. (|2.2|). 
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TIP 



K+ 




u„ 



K+ 



s„ 



B 



)Q 



D 



71^ 



4 




71^ 



Figure 3: The constituent quark loops coupled to K and vr with the insertion of the 
operators Qi or Q2, neglecting soft gluon corrections. (A) gives the leading 0{N^) 
contribution of Qi to the chiral coefficient G^^^, while (B) is the subleading 0{Nc) 
correction to G1^. The opposite happens for Q2 due to the different color structure: 
(C) gives the subleading contribution to G^^ and (D) the leading one to Gl^. 



By computing the two-loop unfactorized diagram in the NDR scheme we obtain 



GM 



I {QQ) j:2 I i-NDR _ Q — 



(5.28) 



By performing a Fierz transformation on the operator 

Q& = -8{sLqR){qRdL) (5.29) 
and using the density building blocks in eqs. (|5.8|)~( |5?TT| ) one obtains, in the NDR 



case, 



M 



In the HV case instead both procedures lead to 



(5.30) 



(5.31) 
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Notice that eq. ( ^.31 ) is not the correct result because of the necessary subtrac- 
tions that must still be implemented in the HV scheme, as discussed in the previous 
subsection. 

Because a similar pattern holds for all ten operators, we have proved that a Fierz 
transformation preserves the xQM result in the HV scheme but not in the NDR one. 
This feature, which was already observed in short-distance calculations [0, 13 1, can 
be understood as a consequence of the prescription of symmetrization of the chiral 
vertices in the HV scheme, which is equivalent to considering the inserted operators 
as four dimensional objects, for which the usual Fierz transformation are allowed. 

As a consequence of these results, when computing matrix elements in the HV 
scheme we will always resort to Fierzing the quark operators in such a way to exploit 
the simpler factorized form. Subtractions are then applied in order to satisfy the 
relevant Ward identities and equations of motion. On the contrary, we do not apply 
any Fierz transformation in the NDR case. We retain only terms up to order M^/A^. 
Ambiguities in the subtraction procedure — like those arising from the expansion of 
the gamma functions in the building blocks of fig. 1(B) — can be shown to be of 
higher order. 



A subtlety arises at this point. The NDR prescription of ref. |11, 13| preserves the 
chiral properties of the operator (Qi — Q2) by means of a special choice of coefficients 
for the evanescent operators. Consistency with such a prescription suggests that we 
have to impose, by an appropriated subtraction, that the operator Qi — Q2 (as well 
as Q9 — Q10) remains a pure octet (8^ x ii?)- As a consequence, the shift in eq. ( |5.25| ) 
is cancelled and the matrix elements induced by Qi and Q2 are the same in the two 
schemes. 



6 The Coefficients of the Weak Chiral Lagrangian 

The results of the previous sections makes it possible to compute the contribution 
of each of the ten operators in eq. ( |2.2| ) to the seven coefficients of the weak chiral 
lagrangian ( f4.2D in both the HV and the NDR schemes. We have included in our 
computation all contributions of order 0{N^), 0{Nc) and 0{asNc). Our result 
depends on the intrinsic xQM parameter Af and three input parameters: /tt, {GG) 
and {qq). 

For the purpose of this computation, it is convenient to rewrite the electroweak 
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operators as follows 



Qi = 


3. ^ 3 
^e,Q, + - 




- ed)AQ7 


Qs = 


3. ^ 3 
^edQe + 2 




- ed)AQ8 


Q9 = 


3. ^ 3 
^edQ3 + 2 


[eu 


- ed)AQQ 


Qio = 


3. ^ 3 
-^edQi + - 




- ed)AQio 



(6.1) 



where 



AQj = 


{sd)Y-A 


(™)v+A 


AQs = 






AQg = 






AQio = 


{S(xdp\^ 


_A (Uf^Ua 



^A • (6-2) 

Prom now on, we identify /^^^ = f = fir- The inclusion in the chiral coefficients of 
the corresponding factors 

- ^VudV*, [ziin) + r2/i(/i)] (6.3) 

is understood. 

6.1 The HV result 

Let us first list the results for the chiral coefficients in the HV scheme. We include 
in the coefficient Gg the contributions of the gluonic penguins, which are pure octets, 
as well as the gluon-penguin-like components of the electroweak penguins (eq. ( |6.lD ). 
We therefore find: 

^8(^3) = /.'^(l-<5(GG> 
Gs{Q4) = ft 
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Gsm = 3e.M/.^(l-6^^ 
GsiQio) = (6.4) 



where 5(gg) is given by eq. (|5.7|) . 

The {V - A)®{V - A) operators Qi,2 and AQg^io yield: 



gIl{Q2) = -ft 



g^l(Qi) = -n 

G\Lm = -j^jt{^-5{GG) 



(GG) 



2' - ^'■'^N^ 
GIl{Qw) = -l{eu-e,)ft 

GldQio) = -l{eu-ed)ft^^{l-6^GG)) ■ (6.5) 

For the constant term which represent the leading contribution of Qs and Qj, 
we find: 

GW(Q7) = -3^(e„-e,)(g<?)2 

GW(Q8) = -3{eu - id) {qqf (6.6) 
while, for their momentum corrections, we have: 

GUm = 3(e„-e,)i-^(g-g) 

GldQs) = 3(e„-ed)^(g-g) 
GiniQy) = -lieu -e,) ft 



GlniQs) = -^(e«-erf)/.'^(l + ^(GG> 



M2 1 f2 
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GiRiQs) = -9^iiu-id)§{QQ) . (6.7) 

6.2 The NDR result 

A similar computation in the NDR scheme leads to a different determination of 
the coefficients because of the shifts discussed in the previous section. In general, a 
shift is present whenever the determination of the coefficient requires the evaluation 

of the unfactorized configuration of fig. 1(A). 
The Gs contributions are now given by: 



X 

M2\ 



X 



3. ,4/. 



For the {V — A) i^i {V — A) operators we have: 
GldQl) = -^/^(l-5{GG)) 

gIl{Q2) = -ft 
gUQi) = -ft 

GlM) = -^/.'(l-5(GG>) 

GIl{Q9) = -lieu-id)ft^{l-S^GG)) 
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GidQ9) = -\{eu-ea)n 



GUQio) = -^(e.-e,)/^^^(l-5(GG>) (6.9) 



For the constant term we find: 



G^'\Qs) = -3{eu - ed) {qqf 1^1 - 3 

G^^HQr) = -3^(^«-^)(^-^)^(l-3|^) (6.10) 
with the corresponding momentum corrections: 

GUm = 3(e.-e,) — ^(g-g)h-3^j 



Glr{Q?>) 


= 3(e„ - 




M2 

V X 


G\r{Q7) 


= -^{eu 






GlniQs) 


= -2(e« 


-e.)/.^^( 




GlniQr) 


X 


(e.-e,)i- 




GIr{Qs) 


X 


(e„-e,)^ 





6.3 Discussion 



(6.11) 



It is perhaps useful to summarize here the approximations we have made in 
obtaining the weak chiral lagrangian of eq. ( |4.2| ) and its coefficients. 

First of all, the lagrangian (|4.2| ) is only given to the O(p^) and any amplitude 
is expected to receive sizable corrections from higher order terms. These have been 
classified [14| but a direct computation of their coefficients, even in such a simple 



model as the xQM, is a rather formidable task. An estimate of the 0{p^) contribu- 
tions to some matrix elements of Qq has been given in ref. ||5|. The NLO corrections 
vary from 10% to 30% depending on M. 

Such an approximation implies that we are also neglecting the effect of the heavier 
meson multiplets. As discussed in section 4, they give in principle rather large 
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contributions in the strong sector, in particular to L5. In the weak sector, their 
effect seems to be less dramatic (see below) — except in the modification they induce 
in the axial coupling between the constituent quarks and the mesons. This is 
clearly a place where future improvement is needed. 

Another approximation we have made is in keeping only the first term in the 
expansion in M^/A^, that is characteristic of the model. M^/A^ is small enough 
to make us confident of our result. Anyway, insofar as the xQM is just a (very) 
simple model, it is not clear whether going to the next order would result in a real 
improvement. 

Finally, the chiral coefficients are given to 0{1/Nc). In this respect it is worth 
noticing that the quark operators Q4 and Qq do not induce any 0{1/Nc) correction 
to the coefficients. This happens because of kinematics and of the flavor singlet 
structure of the currents (which induces cancellations among u and d flavor exchange 
in the subleading configurations). As a consequence, gluonic corrections as well are 
absent for Q4 and Qq, appearing first at 0(1/ N^). Similarly, no gluonic corrections 
appear for the operators Q5 and Qj at 0{1/Nc) because of their color and chiral 
structures. Notice that gluon corrections are in the form (1 — (^(gg>) for LL operators 
and (1 + 6(GG)) ^^i those with LR current structure. 

7 The input parameters and M 

The quark and the gluon condensates are two input parameters of our com- 
putation which are in principle free. Their phenomenological determination is a 
complicated issue (they parameterize the genuine non-perturbative part of the com- 
putation) and the literature offers different evaluations. 

For guidance, we identify the condensates entering our computation with those 
obtained by fitting the experimental data by means of the QCD sum rules (QCD- 
SR) [^3| or lattice computations. In our approach, we vary these input parameters 
within the given ranges in order to obtain a value for the amplitude we are interested 
in. 
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7.1 Gluon Condensate 



For the gluon condensate, the most recent QCD-SR analysis |24], based on e^e 
data, gives the scale independent result 

(^^GG) = (388 ± 10 MeV)^ (7.1) 

TT 

Such a value is consistent with older QCD-SR determinations |25] as well as with 



another recent one that finds [27| 



(^GG) = (376 ± 47 MeV)^ . (7.2) 
vr 

These values are systematically smaller than the central value of the lattice 
result M 



(^GG) = (460 ± 21 MeV)^ (7.3) 
vr 

which however suffers of a systematic error that is difficult to evaluate. We will take 
(17.21) as the range to be explored in numerical estimates. 



7.2 Quark Condensate 

The quark condensate is an important parameter in our computation because it 
controls the size of the penguin contributions, in particular of the leading operator 
Qq. Unfortunately, the uncertainty about its value is large because of the sizable 
discrepancies among different estimates. 

In the QCD-SR approach, the quark condensate can be determined from ^'5(5), 
the two-point function of the hadronic axial current, as [28| 



^5(0) = -(m„ + md){uu + dd) = 2/2^2 (1 _ 5^) 

= (3.2 - 3.3) X 10^ MeV^ . (7.4) 



This estimate agrees with the most recent determination |29] of the parameter 6-,^ 



that quantifies deviations from the PCAC result. Such deviations are of a few 



percents in eq. (7.4) (but larger in the case of the strange-quark condensate). 

These are scale independent values. To obtain the scale-dependent quark con- 
densate, we use the renormalization-group running masses -|- m^, the value of 
which has been estimated at 1 GeV to be [ p9| ] 

rfiu + md = 12 ± 2.5 MeV (7.5) 
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for Ag^o = 300 ± 150 MeV. The error in ( |7.5| ) reflects changes in the spectral 
functions. In our numerical estimates, we will take as input values the running 
masses at 1 GeV given by ([7.5[). Even though our preferred range of Aq^j-, (see 
eq. (^3)) corresponds to Aq^j-, = 400 ± 100 MeV, we feel that we are not making 
too large an error since this determination is not very sensitive to the choice of 
Aqcd- 

By taking the value ( fT.SD , we find for the scale-dependent (and normal-ordered) 
condensate 

{qq) W = - — . X , _ . N , 7-6) 
the numerical values of 

(qq) = -(238 ± 19 MeV)^ (7.7) 
at 1 GeV and 

{qq) = -(222 ± 19 MeV)^ (7.8) 

at 0.8 GeV. The error in eqs. (^)-(|7!|) is due to that in (^. 

On the other hand, a recent determination of the quark condensate in lattice 
simulations with quenched Wilson fermions [^] finds a value of 

{qq) = - (257 ± 27 MeV)^ (7.9) 

at 1 GeV. 



A similar simulation with dynamical staggered fermions |31] yields the rather 
large result 

{qq) = -(380 ± 7 MeV)^ , (7.10) 

which is probably an overestimate to be discarded for our purposes. 

As we can see, the actual for the quark condensate should be varied in the range 

- (200 MeV)3 < {qq) < -(280 MeV)^ (7.11) 

in order to include the central values and the errors of the QCD-SR and lattice 
estimates. 
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7.3 The Constituent Quark Mass M 

Concerning the values of M, we take the point of view that it is an arbitrary 
parameter to be best fitted by comparing the predictions of the model with the 
experiments and by minimizing the theoretical 75-scheme dependence. This point 
will be discussed further at the end of this paper. 

Let us only comment on what values we consider reasonable for M. As a gen- 
eral rule, we expect M to represent the constituent mass for the quarks inside the 
Goldstone bosons (see eq. ( 3.11| )), and therefore to be roughly 



M PS 200 - 250 MeV , (7.12) 

as consistently estimated in processes [^] involving mesons. Such a value is smaller 
than the value M ~ 330 MeV often quoted that originates from baryon physics. 

8 The Matrix Elements 

We wish now to apply our results to kaon decays by computing the amplitudes 
Aoo = a{k^^ ^°7r°) , A+^=A (k° ^ TT+vr-) (8.1) 

and 

A+Q = a(k+^ 7r+7r°) (8.2) 



by means of which is possible to discuss, for instance, the A/ = 1/2 selection rule 
in kaon physics and e' /e. 

8.1 Chiral Loops and Long-Distance Scale Dependence 

The long-distance scale dependence of the matrix elements is contained to this 
order in the one-loop corrections induced by the propagation of the mesons, a con- 
tribution that we have so far ignored. Their effect is important and essential in the 
matching procedure between Wilson coefficients and matrix elements that is central 
in our approach. 

The computation is rather involved because of the many terms; the diagrams 
(see fig. 4) include the vertex as well as the wave-function one-loop renormalization 
(see ref. for the Feynman rules and the details of the computation). 
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We give the final results for the meson loop renormalization of the matrix ele- 
ments for the three processes in eqs. ( ^■lD - (l8.2| ). The contribution to the tree-level 
amplitude of each term of the chiral lagrangian is formally factorized out as a func- 
tion of the input parameters, while the corresponding loop renormalization is given 
in the form of numerical coefficients. These numerical coefficients are complicated 
functions of the masses and the coupling constant. They are made of polynomial 
terms, generally of the order of 

2 

and logarithmic terms of the order of 

2 2 

„,„ ln— |, 8.4 
(47r/)2 ml' 

where the masses can be any among ttt-tt, mx and m,,. The values of the masses and 
other input variables are those given in table 3 of appendix C. 

Let us recall here that the leading-log approximation would be particularly crude 
in this case since the large mass separation between m^ and m,r makes the renor- 
malization group scale particularly uncertain. 

The resulting corrected amplitudes for K and K — > vr^vr are given 

by: 

— nJOA ( n QOA -^ n dddi -a- n ^.^^.^in 
P 

V2, 



aooiQi) = - ^GolQi) (0.904 + 0.444 i + 0.255 In ^- 



{m\ - ml) [ciLiQi) (0.367 + 0.444 i + 0.0747 In 

- GiiiQi) (l + 0.349 0.0184 i 0.135 In 

- GsiQi) (l + 0.716 + 0.463 i + 0.210 In /i^) 
+ GljiiQi) (l + 0.349 + 0.0184 ? + 0.135 In 



and 



a/2 

-Js^l GlniQi) (0.868 + 0.444i - 0.132 In /i^ 

/2 

+j^rnl: GiiiQi) (o.719 + 0.444 i + 0.203 In (8.5) 



/2 

a+-{Qi) = - yjGo(Qi)(l + 0.708 + 0.240 i + 0.203 In 
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Figure 4: One-loop chiral renormalization of the the tttt amphtudes. The 

black box represents the insertion of the weak chiral lagrangian, whereas the black 
circle denotes the insertion of the 0{p^) strong chiral lagrangian. For each chiral 
coefficient more than a hundred diagrams are generated by propagating in all allowed 
ways K, it and 77. 

V2, 



-^{m\ - ml) [GlLiQi) (l + 0.768 + 0.240 i + 0.278 In 1?') 



+ G^LiQi) (0.0525 - 0.222 i + 0.0688 In fi^ 

- GsiQi) (1 + 0.716 + 0.463 i + 0.210 In 

- GiniQi) (0.000127 - 0.222 z + 0.0583 In 

-j^ml GljiiQi) (1 + 4.17 + 0.240 i + 1.29 In /x^^ 
/2 

+j^mj^ GljiiQi) (1 + 3.47 + 0.240 i + 1.40 In f?^ . (8.6) 

where all dimcnsionful parameters must be taken in units of GcV. 

As a useful check we have also computed directly the meson loop renormalization 
for the amplitude — > tt+tt'', which is a pure AI = 3/2 transition: 

a+oiQi) = - -^Go(Qi)(l- 0.196 - 0.204? - 0.05101n /x^) 

-Jsi'^K- ml) [GlLiQi) (1 + 0.402 - 0.204 i + 0.204 In /x^) 
+ G^LLiQi) (1 + 0-402 - 0.204 z + 0.204 In /x^) 
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- GljiiQi) (l + 0.349 - 0.2042 + 0.193 In 
-j^rnl GljiiQi) (l + 3.30 - 0.204 i + 1.42 In fi^) 

+J3^K GiniQi) (l + 2.75 - 0.204 i + 1.20 In fi^) . (8.7) 



The renormalization at any scale // can readily be obtained from eqs. (8^) and 

3- 

As the reader can easily check, our result satisfies all the expected symmetry 
properties, that is 

A+- = Aoo (8.8) 
for the octet amplitudes, and 

^.0 = (8.9) 
for the other parts. 

A similar computation for G^,'^ and Gg is reported in ref. |3^]. Unfortunately, a 
comparison with ref. is difficult because the authors do not report all the details 
of their computation. We find that our imaginary parts are almost identical to 
theirs, while the real parts of the renormalization computed at the scale rrir^, where 
the results of [^] are given, differ, even though they are of the same order. At any 
rate, the fact that eqs. ( |8.8| )-(|8l9D are exactly satisfied for all coefficients makes us 
confident on our results. 

Of course, the polynomial parts of these corrections receive contributions at the 
tree level from the next order terms in the chiral lagrangian. These are controlled by 
the two parameters Ki and K4 defined in reference |34|, where they are estimated 
to be 

Ki = (0.4 ± 1.2) X 10"^^ 

= (0.3 ± 1.4) X W-^^ . (8.10) 

A similar result was obtained for Ki in the framework of the xQM in H]. 

These values correspond to a renormalization of the amplitudes Aq and A2 at 
the percent level, which is an order of magnitude smaller than the renormalization 
induced by chiral loops. Such a small contributions are in agreement with the 
results of the factorization and the vector dominance models, where these coefficients 
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are vanishing. As a consequence, our results are not appreciably modified by the 
inclusion of these terms. 

We recall here that, contrary to the usual treatment, in our approach these 
next-to-leading-order contributions do not cancel the scale dependence of the chiral 
loops. 

On top of this renormalization, we should also include the one-loop determination 
of / in terms of /tt and fx- This is taken into account in our computation by 
replacing / by the one-loop renormalized value 

/i = 0.087 GeV [18] (8.11) 

in the tree level amplitudes. This correction amounts to a further 20 % of renormal- 
ization for the amplitudes. 

Some of the loop renormalizations (the last two terms in eqs. (^)-( |8l7|) ) appear 
to be large when compared to the partial tree level amplitudes that we have fac- 
torized out. This is a notational artifact since they remain always smaller than the 
leading tree level amplitude (~ \/2//^). The overall renormalization is large but still 
under control; as a matter of fact, it is large in the / = channel and small in the 
1 = 2 one (except for the subleading LR momentum corrections), thus distributing 
itself in the right amount to bring the matrix elements of the next subsection closer 
to their experimental values. 

In refs. |9| a similar computation was performed by means of a cut-off regu- 
larization, identifying then the cut-off with the dimensional regularization scale of 
the Wilson coefficients. The two approaches lead to different results. Among else, 
the meson loop renormalization of the amplitude A2 is strikingly different, being 
suppressed in |^ ^ , whereas it is slightly enhanced in our approach. 

8.2 The Matrix Elements 

The matrix elements of all the operators (2.2) can now be computed. By using 



the chiral lagrangian (|4.2D , we can readily generate the contributions with three 
mesonic external states necessary in the matrix elements. To these results we apply 
the one-loop meson renormalizations discussed in the previous section, obtained from 
eqs. (|8.5|)-(|8^) by subtracting the tree-level parts. 

We write directly the expression for the isospin states 

{Qi)o,2 = {27r,I = 0,2\Qi\K^), (8.12) 
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Accordingly, the corresponding one-loop meson corrections are denoted by ao,2(Qi)- 
The Clebsh-Gordan coefficients for the isospin projections can be found in the ap- 
pendix. 

Wc recall that in the HV scheme one may rely on Fierz transformations and 
use the factorized building blocks when computing the chiral coefficients. Fake HV 
anomalies must then be subtracted according to the discussion in section 5. Using 
the HV results of section 6 we obtain: 



{Qi)o 

{Ql)2 
{Q2)0 

{Q2)2 

(Q3)0 

(Q4>0 
(Q5)0 

(Q6)0 
(Q7)0 

{Q7)2 
{Q8)0 

{Q8)2 
{Q9)0 



3 

V2 



-1 + 



X 



V2 



1 



1 



(GG> 



X 



^ (l - S^GG) 



+ oo(<3i) 
+ a2(<9i) 
+ ao(Q2) 

+ a2(Q2) 



^X (l - 6{GG)) + aoiQs) 

X + ao(Q4) 
2 (gg> 



= 2 



X' + ao{Q5) 



2^3 {qq? 

2 im) 



X' + aoiQe) 



J_MLx' 



= Vg 



1 

2 

72 {qq) 



= 2^3 



n Nc 
{qqf 



f3 

J TT 

\ 

2 

T 

f3 



{qq) 
1 



X' 



2Nr 



1 



-X 



^(l + ^(GG>) +02(^8 



+ ao(Q9) 



(8.13) 
(8.14) 
(8.15) 

(8.16) 

(8.17) 
(8.18) 
(8.19) 

(8.20) 

(8.21) 
(8.22) 



(8.23) 
(8.24) 
(8.25) 
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{Q9)2 = 



1 + ^ (1 - S,ac) 



{Q 



10/0 



{Q 



1 

V2 



1 - ^ (l - <J(GG> 



10/2 



1 + ^(1-5(GG) 



+ 02 (Qg) 
+ oo(Qio) 
+ a2(Qio) 



where 



X 



and 



Y = V3U 



ml + 3 m\ 



M2 
A2 

X J 



^.26) 
^.27) 
^.28) 

^.29) 
^.30) 



\gg) is given by (|577| ). 

The one-loop renormahzation of / is taken into account by replacing / with /i 
in the tree-level amplitudes, which amounts to replacing 1//^ with 1//^ multiplied 
by 



1 , o ^'^ 

1 + 3 — ~ 



1.18. 



^.31) 



In the NDR scheme we are not allowed to Fierz transform the quark operators, 
and we must resort to the direct computation of the unfactorized two- loop diagrams. 
On the other hand, we need not worry about chiral anomalies as in the HV scheme 
and no subtraction is required. Using the NDR results of section 6 we thus find: 



(Qi)o 

{Ql)2 

{Q2)0 

{Q2)2 

{Q3)0 
{Q4)0 
{Q5)0 



1 

-X 
3 

3 

3 



1 + ]| (l - hoc) 



id-.. 



Nr 



(GG) 



X 



1 + 



Nr. 



Hgg) 



V2 
3 

^(^'-5(GG>^)+ao(Q3 
X' + ao{Qi) 

2 im) 



+ ao(Qi) 
+ a2(Qi) 
ao(Q2) 
+ a2{Q2) 



X" + ao{Q5 



i.32) 

i.33) 

iM) 

i.35) 

i.36) 
i.37) 
i.38) 
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{Q6)0 
{Q7)0 



{Qsk 



{Q9)0 
{Q9)2 

(Qio)o 

(Ql0)2 



where 



= 2 



{Q7)2 = 



2VS {qq? 



X" + ao(Q6) 



NcMf^ - ■ 2 



Y' + lx + aoiQj) 



iVcM/2 



X" 



V2 {qq) , V2 



2V3 

-2 



(qq)" 

f3 

jqq) 



1-3 



jqq) 



X" 



Y' (l + 5(GG>) + ao(Q8) 



-V24lt y 



2iVe 

{qq)^lJ 

V2 



1 

~2 
V2 



M/2 - 2Af, 
X-^(2X-X' 



X 1 + 



\GG) 



X' 



+ 02(^8) 



X 



1 



+ 02 (Q9) 



2 
V2 



'^X-X'-^{^-^iGG))x 



+ ao(Qio) 



1 + ^ (1 - S^GG)) 



+ a2(Qio) • 



X" = x\^l-9 , Y' = VSU 



+ 3 [niK-m^j—^ 



A2 



(8.39) 



^.40) 



(8.41) 



(8.42) 



(8.43) 
(8.44) 
(8.45) 
(8.46) 
(8.47) 

(8.48) 



{Qi)2 = for z = 3, 4, 5, 6 in both schemes. 

The equations above show the importance of the corrections of 0{asN) (pa- 
rameterized by the value of the gluonic condensate) as well as of the meson-loop 
renormalizations. In the limit S^qq^ — ^ and zero meson-loop renormalization, the 
HV hadronic matrix elements are the same as those found in the l/Nc approach, 
except for the {V — A) iS) {V + A) operators Q5, Qe, Q7 and Qs for which the detailed 
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form of the matrix elements is characteristic of the model employed. For instance, 
by means of eq. ( 3.19| ) we find 



(Q6>o = -4^X (8.49) 
in the \/Nc computation, where we used 



Cl + C2 = - 1 o o (8.50) 



as determined from A{K'^ — > ^'^v^)/A{'K^ Z^^^^/x) 0. Eq. ( ^.49 ) shows the 



quadratic dependence on the quark condensate which is distinctive of such an ap- 
proach. 

To our knowledge the terms proportional to Y and Y' in the matrix elements of 
Qj and Qg have been neglected so far. As an example, in the framework we 

find for the matrix element {Q^)2'- 

(Q,), = V6 + 2^6 Ua, -c,)ml- c^m],] - , (8.51) 



2Nr 



where the absolute values of ci and C2 remain undetermined (only their sum is 
determined by eq. ( |8.50| )). Analogous contributions appear in the matrix elements 
of Qt- The matrix element ( ^.51 ) correponds to the 1/Nc determination of the chiral 



coefficients 

G^iJ = -6%^(ci-C2) and Gi« = -6M!c2, (8.52) 

X X 

to be contrasted to that of the xQM in eqs. ( |6.7D and ( |6.11| ). In particular, the 
xQM determination of and gives 02/ {ci — C2) ~ 0(M^/A^), thus making 
the term proportional to m\ of the same order of magnitude of that proportional 
to m^. 

Although these additional contributions do not affect the estimate of the A/ = 
1/2 rule, which is little sensitive to the electroweak penguins, they do have an impact 



on the determination oie' je [3£] 



8.3 A Final Comment 

In studying any physical process, we must consider the matching of the hadronic 
matrix elements with the Wilson coefficients at a scale typically of the order of A^. 
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At the NLO, the Wilson coefficients of ( p.l| ) depend — beside the energy scale — on 
the 75 scheme employed [11, 13 1. One should thus verify to what extent the physical 
amplitudes turn out to be renormalization scale independent, as a result of a balance 
between short- and long-distance scale dependence. On the other hand, one should 
check the 75-scheme independence of the results. The latter may allow us to restrict 
the range for the constituent quark mass M, as it was already shown in ref. in a 
toy model for e' je. 

These issues deserve a separate and detailed analysis which is presented in 
refs. [^] and |36|. In the first of these two papers we address the computation 
of the A/ =1/2 rule in the CP-conserving kaon decays. The results obtained are 
extremely encouraging, allowing us to give in [^] a new estimate of the CP-violating 
parameter e'/e in the standard model. A similar approach for the AS" = 2 lagrangian 
is also under investigation |37]. 
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A Chiral Quark Model 
A.l Feynman Rules 

The free propagator for the constituent quark is given by 

where ^ = j-p- The same propagator in the external gluon field (fixed-point 



IS 



Slip) = - ^ dsT^'GI^^ (p2 _ M2)2 ' 



where 



(A.l) 
?auge) 

(A.2) 
(A.3) 



and (Tp 

Other useful formulae are: 



Tr glT^T^Gl^G^^p = '—{^GG) i6^aSu/3 - 
where Tr T^T^ = Tr X^X^/A = 5"^ /2, and 
fx'^V^, = 12 I ; cT^"7pC^M- = • 

The relevant meson-quark interactions are derived from the lagrangian 
which we write here as 

M _ M 



(A.4) 



(A.5) 



c 



xQM 



-Mqq + 2i—q-i^Ilq + 2-^qWq 



(A.6) 



where 



n 



and 



2V V2 



^0 


7r+ 


K+ 


n~ 






K~ 







(A.7) 



vr 



V6 



(A. 
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In the case of a single meson interactions one obtains 



^75 Hq = --j= (u'y^u tt^ — dj^dir'^ + ^753 + uj^d tt"*" + • • • j . (A. 9) 

The relevant Feynman rules are therefore given by: 



C/75S = K'^ ■U75S = tt'^ wy^d — coupling: 

7r° c^ygo? — coupling: '^~J'^^ 

TT W75U — coupling: ~~^75 

7r° ds - coupling: -i -^-^ (A. 10) 



TT ds = K uu = tt"*" us — coupling: +i 
K'^ K~ ss = ir~ uu = ir'^ 7r~ dd — coupling: +i 

f 

TT^ US = tt"*" tt" ud — coupling: +i 



M 

? 
M 

? 
M 

M 



7r° 7r° uu = 7r° 7r° — coupling: +i 2 

All meson fields are entering the vertex. The same rules hold for the conjugated 
couplings. 

A. 2 FiERZ Transformations and Clebsh-Gordan Coefficients 

The relevant Fierz transformations (for anti-commuting fields) are the following: 

aa7;i(l± 75)^/3 C/37'*(lT75)c^a = -2 ao,(l T 75)da C/3(l ± 75)6^5 
aa7/i(l±75)^'/3C/37'*(l±75)da = aQ7ju(l ± 75)^0 c/37''(l ± 75)6/3 . (A.ll) 

We also have for SU{Nc) 

6ap6^5 = -^S^sS^f, + 2 T^sT^^ . (A.12) 



c 



The SU{2) Clebsh-Gordan projections are as given by 

i(^+_-^oo) =/|^+o. (A.13) 



A 
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The SU{3) projections are [39| 

|8, i) = Tr (AaS^'D^s) Tr (AfS^'D^S 



- Tr (A^StD^Sj TV (^\\^^D^'^j 
|27, i) = Tr (A^StLl^s) Tr (AfS+D^s) 

+ 4 Tr (a^S+D^s) Tr (aJS+D'^s) 
+ 5 Tr (A^St^^s) Tr (A^StZJ/'E) 
|27, I) = Tr (A^S"fL>^s) Tr (AfS+D^s) 

+ Tr (A^StD^s) TV (aJS+D^s) 

- Tr (A^StD^s) TV (A^StL>^'s) . (A.14) 
Therefore, we have 



|27) = ^27, i) + - |27, i) 



2 

= - Tr (a^S+I^^s) Tr (a^S'^Z^'^E 

+ TV (AiS+D^s) Tr (aIs+D^'s) . (A.15) 

A. 3 Dimensional regularization 

We work in the MS scheme. In the naive dimensional regularization (NDR) 
everything is continued to d dimensions and the same four-dimensional rules applied. 
We therefore have that 

{7^,74 = 25^, (A. 16) 

and 

< = 'i- (A.17) 

The 75 matrix is defined so as to anti-commute in any dimensions as 

{7m.75} = 0. (A.18) 

The term naive refers to the fact that such a prescription leads to manifest 
algebraic inconsistencies. 
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In the 't Hooft-Veltman rcgularization (HV) the Dirac matrices are separately 
considered in 4 (tilded quantities) and d — 4 (hatted quantities) dimensions so that 



(A.19) 



The two sub-spaces are orthogonal to each other: 



(A.20) 



and 



9i: = d-A g^ = 4 and « = 0. 



(A.21) 



The 75 matrix is defined as anti-commuting in 4 dimensions and commuting in 
d — 4; therefore 



In both schemes the external momenta are kept in four dimensions. Chiral cur- 
rents must be symmetrized in order to have a unique definition. This is immaterial 
in the NDR case but gives 



in the HV case. 

A. 4 Trace Formulae 

In writing the chiral AS = 1 lagrangian we have rewritten single traces as the 
product of two traces. The trace factorization properties can be easily shown. The 
eight SU (3) Gell-Mann matrices A" together with the identity matrix A'^ form a basis 
for the 3x3 complex matrices with non- vanishing trace. The standard normalization 
Tr A'^A^ = 2(5"^ implies 



{7m'75} = [7^, 75] = 0. 



(A.22) 



The rules above lead to 



{7^, 75} = 2 75r- 



(A.23) 



+ 75)7^(1 - 75) = %(1 - 75) 



(A.24) 




(A.25) 
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Given two complex 3x3 matrices A and B we can write 

8 



a=0 



X^B = ^A^6^ (A.26) 



a=0 

where a'^ and are complex numbers. From eqs. ( |A.25D -( |A.26 ) it follows that 

8 

2 



1 * 

Tr (^X^A X'^b) = 2 II Tr (a^'AJA) Tr (A-^A^S) . (A.27) 

a=0 

Finally, by using the explicit form of the Gell-Mann matrices and the identity, we 
have, for the cases of interest, 

TV (xlA Ajs) = TV (X\A^ Tr (a?^) (A.28) 

Tr (X\A X\B^ = Tr (x\A^ Tr (a^^) . (A.29) 
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B Chiral Perturbation Theory: Feynman Rules 



We give in the two following tables the chiral perturbation theory rules needed 
in computing the coefficients of the weak chiral lagrangian. Many more rules are 
necessary in the chiral loop computation (see p3[). 





Feynman Rule 


Coefficient 


i^°(pi)vr°(p2) 









if. 


Gs 


iyl Pi ■ P2 


-ijiPi ■P2 







iji Pi ■ P2 


r^b 


ijT Pi ■ P2 





Glr 





-ijiPi ■P2 


'^LR 


-ijipi ■P2 





Glr 











Feynman Rule 


Coefhcient 


K\vi)^\P2y\P3) 


K0(pi)^+(P2)^-(P3) 


Gs 

ria 
^LL 

rib 
^LL 

ria 
'^LR 

rib 
'^LR 

ric 
'^LR 




Jlp [2^51 • (P2 +P3) P2-P3\ 



Jlp [2^*1 ■ (P2 +P3) P2 -PsJ 



Jlp [2PI • +P3) P2 -^3] 




-4 

^P2 ■ iP3 -Pl) 
-7X^3 • iP2-Pl) 

-ylpi • {P2-P3) 

JtP3 ■ {Pl +P2) 

-^rPi ■ {P2-P3) 

-'^Pl-P2 



Table 1: Relevant Feynman rules of chiral perturbation theory. All momenta are 
entering the vertex. 
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C Input Parameters 



parameter 


value 


/tt = fn+ 


92.4 MeV 


Ik = fK+ 


113 MeV 




138 MeV 


rriK = rrixo 


498 MeV 


rur, 


548 MeV 




350 ± 100 MeV 


rnu + rrid (1 GeV) 


12 ± 2.5 MeV 




-(200 - 280 MeV)^ 


{asGG/n) 


(376 ± 47 MeV)^ 



Table 2: Table of the numerical values used for the input parameters. 
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